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Recent precision measurements of mesoscopic persistent currents in normal-metal rings rely on 
the interaction between the magnetic moment generated by the current and a large applied mag- 
netic field. Motivated by this technique, we extend the theory of mesoscopic persistent currents to 
include the effect of the finite thickness of the ring and the resulting penetration of the large mag- 
netic field. We discuss both the sample-specific typical current and the ensemble-averaged current 
which is dominated by the effects of electron-electron interactions. We find that the magnetic field 
strongly suppresses the interaction-induced persistent current and so provides direct access to the 
independent-electron contribution. Moreover, the technique allows for measurements of the entire 
distribution function of the persistent current. We also discuss the consequences of the Zeeman split- 
ting and spin-orbit scattering, and include a detailed and quantitative comparison of our theoretical 
results to experimental data. 



PACS numbers: 73.23.Ra, 73.23.-b, 05.30. Fk 

I. INTRODUCTION 

Persistent currents in normal-metal rings threaded 
by an Aharonov-Bohm flux constitute a paradigm of 
quantum-coherence effects in the thermodynamic prop- 
erties of mesoscopic systems. While the history of per- 
sistent currents dates back to the early days of quantum 
mechanics^ and of superconductivity,^ they were studied 
intensively starting with the seminal paper by Buttiker, 
Imry, and LandauerP 

Most experiments to date detected persistent currents 
using SQUIDs (superconducting quantum interference 
device) as magnetometers. A different technique was re- 
cently developed by Bleszynski-Jayich et al™ which is 
much more sensitive and allows for precision measure- 
ments with lower back action and over a wider range 
of magnetic fields. The high-precision cantilever torque 
magnetometer relies on the interaction of the magnetic 
moment associated with the persistent current and a 
large applied magnetic field. This interaction shifts the 
resonance frequency of a microcantilever on which the 
rings are located. Measurements of the frequency shift 
allow one to extract the persistent current quantitatively. 

This paper extends the existing theory of mesoscopic 
persistent currents to include a large applied magnetic 
field. We focus on metallic samples with diffusive elec- 
tron dynamics for which the applied magnetic fields are 
non-quantizing. Our results hold for both normal-metal 
rings as well as rings made of nominally superconducting 
materials (provided that the magnetic field significantly 
exceeds the superconducting critical field H C 2) and in- 
clude the effects of spin, namely Zeeman splitting and 
spin-orbit scattering. 

Within an independent-electron model, the flux- 
periodic persistent current is strongly sample specific 



with both magnitude and sign depending on the details 
of the disorder configuration and the geometry of the 
ring. As a result, its ensemble average (I) is, even in the 
canonical ensemble small compared to its second mo- 
ment (J 2 )02lso that the latter describes the typical mag- 
nitude of the persistent current of an individual ring. The 
typical persistent current is 0o _ periodic and, in a diffu- 
sive metallic ring, has an amplitude of the order of e/rrj, 
where td denotes the diffusion time of an electron around 
the ring. 

The disorder-averaged persistent current is dominated 
by the contribution of electron-electron interactions^ 
and is </>o/2-periodic. In a normal-metal ring, it is of the 
order of A(e/T£>), where A is an effective electron-electron 
coupling constant. While A is of order unity in lowest or- 
der perturbation theory, higher-order contribut ions are 
expected to reduce its magnitude significantly!^^ In 
rings made of superconducting materials, a related mech- 
anism leads to a current due to superco nduct ing fluctu- 
ations above the critical temperature T e ! 13 * 14 l 

These theoretical expectations have been test ed 
in several experiments, including metallic] 4 * 16 * ^ 
semiconducting,' 2 ^ as well as superconducting^ rings. 
While results of early experiments with metallic rings 
were in apparent strong disagreement with theoretical 
predictions, a more recent SQUID-based experiment!^ 
yielded data reasonably close to theory. Finally, the 
measurement of the typical persistent current reported 
in Ref. El agrees, without any adjustable parameters 
and over a wide range of experimental variables, with 
the predictions of the model of nonintcracting diffusive 
electrons described here. We also note in passing that 
there is a closely related set of works, both experimental 
and theoretical, which explores the magnetic response 
of singly-connected mesoscopic systems, see e.g., Refs. 
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Persistent currents have also motivated a multitude of 
further theoretical considerations. Among other results, 
it was suggested that the persistent current is highly sen- 
sitive to a variety of subtle effects, including the cou- 
pling of the ring to its electromagnetic enviro nment 25 * 26 ! 
as well as magnetic impurities within the ringJ^EU This 
indicates that accurate measurements and understanding 
of persistent currents in various settings would address a 
number of interesting questions in many-body condensed 
matter physics. 

This paper is organized as follows. In Sec. [TTJ we dis- 
cuss the flux dependence of the persistent current. Sec. 
|III| discusses the effects of the strong magnetic field on 
the persistent current within the independent-electron 
model, including the effects of the Zeeman energy and 
spin-orbit scattering. Sec. |IV| focuses on the interaction 
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contribution to the persistent current. Sec. |V| contains a 
detailed comparison between our theoretical results and 



In this case, gauge invariance implies flux periodicity, 



1(4) = 1(4 + fa), (I) 



where the period is given by the flux quantum c/)a = h/e, 
and time-reversal invariance gives the relation 



(2) 



the experimental data of Ref. [31 We conclude in Sec. VI 



As a result, the persistent current vanishes at integer and 
half-integer multiples of the flux quantum and can be ex- 
pressed as a Fourier series I(4>) = YlpLi Ip sin(2np4/ 4o) ■ 

It is instructive to deduce the consequences of this 
Fourier decomposition for the current-current correlation 
function 



II. FLUX PERIODICITY OF THE PERSISTENT 
CURRENT 

Conventionally, persistent currents are discussed in the 
limit of a pure Aharonov-Bohm flux threading the ring. 



C(<f>,<p') = (1(4)1(4'))- 



(3) 



Here, (. . .) denotes a disorder average. We anticipate 
that the Fourier components I p are mutually uncorre- 
cted, i.e., (Iplp') = (Ip}5 pp >. Then, we obtain 



oo oo /r2\ 

C(0>') = ^(7 p 2 ) S in(2 7 r^/0o)sm(27r^7^ o ) = X] 1 f Z 



P =i 



{cos[2^p(0 - - cos[2 7 rp( ( /» + </>')/ <h]} ■ (4) 



P =i 



r 



Within the diagrammatic approach to diffusive electronic 
systems, the two terms depending on (4 — cfj') and (<fj + 
</>') have immediate interpretations as the diffuson and 
cooperon contributions, respectively.^ Both contributions 
are of the same magnitude but depend differently on the 
magnetic flux. 

In the presence of an additional large magnetic field B 
penetrating the metal ring, one expects that the cooperon 
contribution is strongly suppressed. This leads to a 
change in the flux dependence of C(4, 4') which can 
also be obtained directly from symmetry considerations. 
While gauge invariance and hence the flux periodicity 
persist, the additional magnetic field changes the time- 
reversal relation into 7(7?, 4) = —I(—B, —<f>). As a result, 
the current is no longer odd in the Aharonov-Bohm flux 



4 alone, and the Fourier series takes the more general 
form (at fixed B) 

oo 

= cos(2ttp</)/0o)+7(-} sin(2^/0 o )}. (5) 

p=i 

If we again anticipate that the Fourier components are 
mutually uncorrelated, 



(i^i^) = o 



(G) 
(7) 



and that, moreover, ([7^] 2 ) = ([7^ ^] 2 ), we find 



C(<M') = ^([7W] 2 )[sin(2^/0 o )sin(27r P( /)70 o ) +cos(27T^/0 o )cos(2^70 o )] 
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^([/W] 2 )cos[2^(0-0')/0o]. 
p=i 



(8) 




For non-interacting electrons, the (grand-canonical) ther- 
modynamic potential 51 can be expressed as 



FIG. 1: (a) Diffuson and (b) cooperon diagrams for the au- 
tocorrelation function of the density of states. Full lines rep- 
resent electronic Green functions, dashed lines correspond to 
disorder scattering. 



In agreement with expectations, our analysis implies that 
in the presence of a large magnetic field B, the current- 
current correlation function has the flux dependence of 
a diffuson contribution. Note that the magnitude of the 
persistent current, (I 2 (</>)), becomes independent of flux. 
As a special case, this also implies that the persistent 
current can be nonzero at zero flux. 

It is interesting to note that in the presence of a large 
magnetic field, the flux dependence of the persistent cur- 
rent can also be written as 



B) = —T J dE v{E, B) ln[l + e -K B -fi) (11) 

in terms of the density of states v{E,H). Here, f3 — 1/T 
denotes the inverse temperature. (We use units ks = 1 
and h = 1.) Here, the magnetic field B includes both 
the Aharonov-Bohm flux 4> threading the ring and the 
magnetic field penetrating the ring. For definiteness, we 
will from now on decompose the full magnetic field into a 
pure Aharonov-Bohm contribution and an in-plane field 
.By penetrating the ring. Accordingly, we will drop the 
vector nature of B in the following although it should be 
kept in mind that in principle, the persistent current is 
not an isotropic function of magnetic field. 



The thermodynamic potential f2(/i, B) a ^ finite tem- 
perature can be related to its zero-temperature limit 



!(<P) =^2 t p cos(2npcl)/<f> - a). 



(9) 



Comparing with Eq. rt5 



I p cos a and 



r(+) 



yields the identities Ip 
I p ' = IpSiria. Then, we automatically re- 
produce Eqs. ([6| and (7j) by assuming that the phase 
offset a has a uniform distribution over the disorder en- 
semble. This also yields the relation (1%) = 2([/p ± ' ) ] 2 ). 

In the next section, we verify these flux dependencies 
explicitly within the model of diffusive non-interacting 
electrons. 



Q (fi,B) 



dE(E — n)v(E, B) 



(12) 



as 



III. INDEPENDENT-ELECTRON 
CONTRIBUTION 

Current-current correlation function 



SI(H,B) 



dE 



dE 



n (E,B) 



(13) 



The persistent current is obtained as the flux- 
derivative of the thermodynamic potential 



I = — 



on 



(10) 



in terms of the Fermi-Dirac distribution fJE). Thus, 
the current-current correlation function Ci{B 1 B') = 
(I(B)I(B')) takes the form 



d(B,B') = J dEdE' 



dU(E)\ ( dfp(E') 



dE 



dE' 



(n (E,B)n (E',B / )) 
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^^( l-expt-w )^'^ 5 ^' 5 '' (14) 

Here, we used in the second identity that the correlator depends only on the energy difference e = E — E' so that we 
can perform the integral over the sum a = E + E 1 . Thus, we are led to consider the zero-temperature autocorrelation 
function Cj(E, B; E', B') = (I(E, B)I(E', B 1 )) of currents at different chemical potentials E and E' as well as fields 
B and B' . 

Within a model of non-interacting, diffusive electrons, the calculation of 

cf ) {E,B-E',B')= ( dE 1 j dE 2 (E 1 -E)(E 2 -E') 1 ^—(v(E 1 ,B)i;(E 2 ,B')) (15) 

starts from the familiar diagrams in Fig. for the disorder-averaged autocorrelation function of the density of states.^ 
Note that both the diffuson and the cooperon diagram contribute to the persistent current. The diffuson diagram 
depends on the difference A_ = A — A' of the magnetic vector potentials, the cooperon diagram on the sum 
A + = A + A'. Performing the integrations over the fast Green-function arguments in the diagrams of Fig. [IJ one 
arrives at the expression 



C?\E.B-E',B') = -1-VReV [ 



/>oo />2cr 






da de 


M. 


Tr ^ 


Jo J -2a 





, -D[V - ieA ± } 2 +i{e + E - E 1 ) 



(16) 



Rewriting the square of the diffusion pole as a derivative with respect to e and integrating by parts yields 

1 <9 2 f°° f 2a f 1 \ 

cf\E,B ] E\B') = -^——Y dalm rfeeTr — -= — — - 9 — . (17) 

7 4tt 2 <90<90' ^W ]_2a \-D[V -ieA ± } 2 + i{e + E - E') J 



Here, D denotes the diffusion constant and we limit at- 
tention to spinless systems. (Effects of spin will be dis- 
cussed separately in Sec. Ill C ) 

In Eq. (17), the trace is over a space of wavefunctions 



-0 satisfying the condition 



n • 



S A±]V>| S - 0. 



(18) 



at the surface S of the metallic ring, (n denotes denotes 
a unit vector normal to the surface.) In general, this 
boundary condition makes the evaluation of Eq. (17 1 a 
tedious problem. 

To simplify this problem, we use a model in which 
the in-plane magnetic field is taken to be of constant 
magnitude and to point along the azimuthal direction 
around the ring. While this toroidal-field model is clearly 
different from experimental realizations, we expect that 
it gives a qualitatively and, for certain quantities, even 
quantitatively correct account of the consequences of a 
large magnetic field penetrating the ring. Specifically, 
we expect that the predictions for the correlation field B c 
are parametrically correct while the numerical prefactor 
would reflect the particular field configuration. At the 
same time, predictions for the typical current amplitude 
will bc quantitatively correct because the large in-plane 
field drops out of the final expressions. 

Some considerations for more general field configura- 
tions are collected in an Appendix. 



B. Toroidal magnetic field 

The simplification of the toroidal-field model derives 
from the fact that in this case, the eigenvalue problem 



D[V - ieA±] 2 V> = Sip 



(19) 



together with the boundary condition in Eq. (|18|) can be 
solved by separation of variables. Let us consider a ring 
defined as a cylinder of length L (along the ^-direction) 
and radius R (in the x — y-plane) with periodic boundary 
conditions in the z-direction. The total vector potential 
A is a sum of the Aharonov-Bohm contribution Aj^ = 
(<p/L)z describing the flux threading the ring and the vec- 
tor potential A|| = (B\\/2)z x r of the in-plane magnetic 
field penetrating the ring. Then, the eigenvalue problem 
in Eq. (19) separates with ip(x,y,z) = x( x i u) exp(ifcz) 
where 



S = E c (n - ip±) 2 + e± 
with n = 0, ±1, ±2 . . . and 

- DM. - "fyf + (d y + l ^ x f] x 
Here, we defined the Thouless energy 



e_LX- 



L 2 



(20) 



(21) 



(22) 
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and the dimensionless flux variable <p± = <fr±/4>o- Note 
that in order not to introduce unnecessary numerical 
prefactors into equations, this definition of the Thouless 



energy differs by a factor of four from the definitions em- 
ployed in Refs. H]and[9l 



Inserting these eigenvalues into Eq. (17), we find 



± cj_ n u 



2cr | 



(23) 



Performing the sum over n by Poisson summation and measuring all energy variables in units of the Thouless energy, 
one obtains 



cf\E,B;E',B' 



2tt 



F 2 pp. 00 /-oo f 2a exp(-27rp A /e| + i(e + E - E')) 

7EEE cos ( 2 ^±)/ d ^ d ^ i ( 24 ) 

' ± ex p=i Jo J - 2rT Je±+i{e + E-E') 



The integrals over e and a can be readily done to yield 

,(0),- „ SE 2 d 2 



Cf'(E,B;E',B') 



7 E E E c °s( 27 rpip ± )F p (z ± ) 



(25) 



where z± = [ej_ + i(E — E')]/E c and where we defined the function 



F p (z) = Re 



(2ttp) 5 (27rp) 4 (2ttp) 3 
I 



(26) 



We note that this result is valid for spinless fermions. 
Effects of spin will be discussed below in Sec. HIC 
It is interesting to compare the result in Eq. (25) with 



the corresponding correlation fun ction for the conduc- 
tance fluctuations of a metallic ring! 31 * 32 ^ Indeed, the flux- 
sensitive contributions to the correlation function of the 
conductance at different magnetic fields differ from our 
result for the persistent current (apart from an overall 
prefactor) only by the preexponential factor in the func- 
tion F p (z). 

In the absence of the in-plane magnetic field, we need 
to retain only the lowest transverse eigenvalue = to 
exponential accuracy in 2L/R. Then, we find 



(WW)) = 



6Ei 



OO 

E 

p=l 



sin(27rp<p) sin(27rpy/) (27) 



for the current-current correlation, which reproduces the 
result obtained in Ref. |U 

In the limit of a large in-plane magnetic field, the 
cooperon contribution is strongly suppressed since time 
reversal symmetry is broken. This can be seen explicitly 
by computing the lowest transverse eigenvalue e ± pertur- 
batively in 73, for both the cooperon and the diffuson con- 
tributions. This perturbative approach is valid as long as 
R <C &Bi where Ib has to be evaluated for the appropri- 
ate in-plane magnetic fields entering the cooperon (+) 
and diffuson (— ) contributions. (Here, £g = (I/g^h) 1 ^ 2 
denotes the magnetic length.) Due to the boundary con- 



dition of zero normal current, the ground state wavefunc- 
tion |gs) of Eq. (|21|) at zero Bn is a constant with zero 
transverse eigenvalue. Thus, the leading correction to the 
eigenvalue is given by 



ex 



gs 



D 



\De 2 {B\ ] 



-(x 2 +y 2 ) 



and we find that 



327T 2 



LR 

l B 



(28) 



(29) 



For the cooperon contribution, the magnetic field is of the 
order of twice the applied magnetic field. Thus, by Eq. 
( 25 1 this contribution is exponentially suppressed once 



the relevant in-plane field is larger than one flux quantum 
penetrating the ring. 

We first focus on the typical persistent current at zero 
temperature. In this case, the effective in-plane field 
vanishes for the diffuson contribution, while it strongly 
suppresses the cooperon contribution. Thus, assum- 
ing from now on that B\\ is sufficiently large to make 
e_i_(273) ^S> E c , we need to retain only the diffuson con- 
tribution and obtain 



OO _^ 

cos(2Trp[ip - ip']). 



(30) 
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FIG. 2: (Color online) Current-current correlation function FI G. 3: (Color online) Current-current correlation function 

(/((/>, B||)7(0,B|| + AS||)) (in units of (E c /(f>o) 2 ) vs. ASy at 
finite temperatures. The curves are normalized to their value 
at AB|| = and correspond to T = 0.01, 0.02, 0.1, 0.2, 0.5 x E c 
(from bottom to top). The inset shows the same curves but 
plotted logarithmically along the vertical axis. 



perature (solid line). The dashed and dotted lines correspond 
to the contributions from the first and the second harmonics, 
respectively. The inset shows the same curves but plotted 
logarithmically along the vertical axis. 



Comparing with Eqs. ^ and Q, we find 



(+)i 2 \ 



)12\ 



a, 



Ml 



(31) The correlation function can be readily evaluated by com- 



for the harmonics of the persistent current. 

Equation (25), (26 1, and (29) also imply that the cor- 



relation function of the persistent current at different val- 
ues of the in-plane magnetic fields falls off exponentially 
with the magnetic-field difference once the in-plane field 
changes by more than a flux quantum through the cross 
section of the ring, i.e., on the scale of the correlation 
field 



7T LK 



(32) 



Note that the functional dependence of the correlation 
field on L and R remains the same for much more gen- 
eral field configurations but that the numerical prefactor 
in Eq. (32) is specific to the toroidal-field model. A plot 
of the correlation function (/(</>, Bn )I(4>, -By + AB»)) is 
shown in Fig. [2] Its exponential fall-off has important 
ramifications in experiment. The decay of the correla- 
tion function implies that measurements of the persistent 
current at in-plane fields which are significantly separated 
from each other on the scale set by B c are statistically 
independent. We are thus led to the ergodic hypothesis 
that averaging over a sufficiently wide range of in-plane 
fields is equivalent to averaging over the disorder ensem- 
ble. This observation is particularly pertinent in view of 
the novel technique of measuring persistent currents em- 
ployed in Ref. 0] which allows one to obtain the persistent 
current over a wide range of in-plane magnetic fields. 

We close this section by discussing the temperature de- 
pendence of the persistent current at large in-plane mag- 
netic fields. At finite temperatures, the persistent current 
correlation function depends on AB and temperature T 
via the two dimensionless variables, AB/B C and T/E c . 



bining Eq. (14) with Eq. (|25|). Performing the remaining 



integral numerically, we obtain the results shown in Fig. 
[3] for the current-current correlation function and in Fig. 
[4] for the temperature dependence of the typical current. 
We see from Fig. [4] that the temperature dependence can 
be approximated as exponential with reasonable (though 
uncontrolled) accuracy. (Numerical values of the fit are 
quoted in the figure caption.) Moreover, we observe that 
the typical persistent current becomes rapidly dominated 
by the first harmonic as temperature increases. 
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0.00 



T/Ec 

FIG. 4: (Color online) Temperature dependence of the typical 
current (I 2 (</>)) (blue line). The dependence can be well fitted 
by an exponential {I 2 } « c{E c /(f>o) 2 exp(— aT/E c ) with c = 
0.036 and a = 8.2 as shown by the red (dashed) line. The 
dotted and dash-dotted lines correspond to the contributions 
from the first and the second harmonics, respectively. The 
inset shows the same curves but plotted logarithmically along 
the vertical axis. 
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C. Effects of spin 

In weak magnetic field and in the absence of spin- 
orbit scattering, spin enters the persistent current simply 
through a degeneracy factor of two. Thus, Eq. (27) is 
modified into 



2\Ei 



\ " 



— sm(2Tr p(p) sin(2Tr pip'). (33) 



This result includes both the diffuson and cooperon con- 
tributions. 

In a large applied magnetic field, but still without spin- 
orbit scattering, the cooperon contribution is suppressed 
and we have to take the Zeeman energy into account. 
The corresponding spinless result was given in Eq. (30). 



We can include the spin and Zeeman energies by writing 
the persistent current as a sum of the contributions of 
spin-up and spin-down electrons, 1 = 1-^ + 1^. Once the 
Zeeman energy becomes large compared to the Thouless 
energy, there are no correlations between 7+ and Tj, and 
as a result, we find 



^4>l ±J p 



CO _^ 

cos(2Trp[tp - ip']). 



(34) 



The recent precision measurements^ of the persistent 
current were performed on samples whose spin-orbit scat- 
tering length is smaller than or of order of the circum- 
ference of the rings, as deduced from weak-localization 
measurements. For this reason, we now turn to a more 
thorough discussion of the consequences of the electron 
spin, which in addition accounts for the spin-orbit scat- 
tering. This can be done by a standard extension of the 
diagrammatic technique for diffusive systemsPS To be 
specific, we focus on sufficiently large magnetic fields that 
the cooperon no longer contributes significantly. Exten- 
sions to include the cooperon contribution at weak fields 
would pose no additional complications. 



Including spin indices, we define the diffuson 
V s } s f (r, r', e) as shown in Fig. [5] and view it as a 4 x 4- 

S 1 S 2 I I 

matrix D(r,r',e) where (si,s\) labels the rows and 
(s2,s' 2 ) the columns. With the ordering (s,s') = (tt 
, tl- It- II) , onc obtains the equation^ 



-D (V - ieA. 



it 



Hz + iJ s , 



D(r,r',e) 



2tt7V(0)7 



? 5(r- 



(35) 



by the standard procedure, starting with the diagrammatic representation shown in Fig. [5} (N(0) denotes the density 
of states at the Fermi energy and r is the elastic scattering time.) Here, the contribution of the Zeeman energy Ez 
yields the term 



H z = 



while spin-orbit scattering is included through 



H Sl 



2 

3t s . 





-2iE z 

2iE z 





(36) 



1 
2 




o -: 
o o 

2 



-10 1 



(37) 



in terms of the spin-orbit scattering time r so . 

By retracing the steps leading up to Eq. (17) in the presence of spin effects, we obtain for the correlation function 
of the persistent current, 



Cf\E,B;E',B') 



1 



d 2 



i 



4^ 2 d<t>d4>' Jo d<7lm J- 2a de ^ \-D[V - zeA_] 2 + i(e + E - E>) 



Hz + H M 



(38) 



where Tr now denotes a trace over configuration space and the four-dimensional spin space. 

In the limit of large Zeeman splitting, Ez ^S> E c , the modes tt and |t are exponentially suppressed. For negligible 
spin-orbit scattering, we then obtain two massless modes tt ± 11- As a result, the correlation function is twice larger 
than the result for spinless electrons given in Eq. (25), in agreement with Eq. (34 1. As the spin-orbit scattering 
increases, only the density mode tt + tt remains massless and in the limit of strong spin-orbit scattering, we recover 



the result in Eq. (25) for spinless electrons. 



More generally, we can discuss the crossover between the limits of weak and strong spin-orbit scattering rate. One 
finds 



cf\E,B;E\B') = - 1 ° 



At: 2 



dalm 



2n 



de < 



2a 



xTr 



-D[V - ieA_] 2 + i(e + E - E') -D[V - ieA_] 2 + i(e + E - E') 



3r S( 



+ 



4 
3t sc 



-£>[V - ieA_] 2 + j(e + E - E' + 2E Z ) + ^- -D[V - zeA_] 2 + i(e + E - E' - 2E Z ) 
where the trace is now over configuration space only. Specifying again to the toroidal-field model, we obtain 

cf\E,B;E',B') 



8E 2 8 2 ^ 

- 2^cos(27rpv5 ± ) 



7r d<j)d<j)' 



, E-E' eT 



+F p i 



E-E' e ± 
E r + ~ 



4 
3t so 



E-E' + 2E Z 

E r 



E-E 1 - 2E Z 
E r 



C -L 3^ 



(39) 



(40) 



where the function F p (z) had been defined in Eq. (26) 



Combining Eq. (40) with Eq. (14) and setting = 0, 
we can obtain the crossover of the typical current be- 
tween the limits of weak and strong spin-orbit scatter- 
ing for arbitrary temperature. (Note that the results for 
the typical current are not restricted to the toroidal-field 
model.) Corresponding numerical results in the limit of 
large Zeeman splitting (where the last two terms in the 
square bracket in Eq. | |40| ) can be neglected) are plotted 
in Fig. [6j which show that the crossover becomes slower 
as temperature increases. 



IV. INTERACTION CONTRIBUTION 

We now turn to a discussion of the interaction- 
contribution to the persistent current in high magnetic 
fields. Adapting the first-order correction in the interac- 
tion V derived in Ref . [10] to the case of a finite magnetic 
field, one finds for the disorder-averaged contribution to 
the grand canonical potential 



Afl = 



N(0)V 



1 E 
dE coth( — )E 



Si 




*2 Si . 


*2 




1 — ► — 






V 


D 


r' = 




s[ 


* 


s' 2 s'i i 


S'l 



Si 



+ r 



s' 



S'l 



FIG. 5: Diagrammatic representation of the equation of mo- 
tion for the diffuson T> s } a J (r, r', e). Full lines represent elec- 

S 1 S 2 

tronic Green functions and dashed lines denote disorder and 
spin-orbit scattering. 



xReTr- 



•D(V - 2ieA) 2 +iE 



(41) 



Here, V is the Fourier component of the screened 
Coulomb interaction potential averaged in momentum 
spaced In a field much stronger than the upper critical 
field of the ring we may constrain considerations to the 
lowest-order correction, Eq. (41). To estimate the inter- 



action contribution to the average persistent current, we 
again employ the toroidal-field model introduced in Sec. 
Ill Then, the eigenvalue problem and boundary condi- 
tions for the cooperon here are identical to those in Eqs. 



(18) and (19), respectively. 



We denote the cooperon eigenvalues by t^ m ± with 
I, n, m being the radial, longitudinal and azimuthal quan- 




FIG. 6: (Color online) Crossover of the typical current {I 2 ((f>)) 
as function of the spin-orbit scattering rate. The curves, cor- 
responding to temperatures T = 0.01, 0.1, 0.3, 1.0 x E c (from 
bottom to top), are normalized to the value of {I 2 {4>)) in the 
limit of vanishing spin-orbit scattering rate. All curves are 
plotted in the limit of large in-plane field where the cooperon 
contribution is suppressed and the Zeeman energy is large 
compared to the Thouless energy. 
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turn numbers, respectively. Due to cylindrical symmetry 
the cooperon modes can be found by separation of vari- 
ables, with the replacement n — > n — 2(f>/4>o added to 
take into account the Aharonov-Bohm flux. In distinc- 
tion from Sec. Ill, the vector potential A in Eq. (41) 
corresponds to the total field so that Ib -C R- In this 
limit, the radial equation can be approximated to lowest 
order in IsjR as 

d 2 



D 

' B 



dx- 



+ x) 2 + ^gfc„,. 



(42) 

where x — r /Ib is a scaled distance from the center of the 
cross section, n m = mis/R and k n ^ = 2Tr(n—2(j)/(f)o)/L. 
Note that the ratio between the radial and the longitudi- 
nal terms in Eq. (42) is dominated by L/ts- The eigen- 



values can be written as 

2 



n,m,<p 



( L 



A/(K m ) 



where the values of Ao(kv. 



o /2 y \2tt£ b , 

(43) 

for the lowest branch of eigen- 
states (/ = 0) can be estimated by using the variational 
method with a Gaussian trial solution. The function 
A (k) has a shallow minimum Aq = (1 — 2/tt) 1 ^ 2 at n* m — 
(tt 2 — 27r)~ 1 / 4 . Using the eigenvalues , to evaluate 



the trace in Eq. (41), it is straightforward to show that 
the contribution to persistent current AI — —dAQ/d(j> is 
periodic in <fi —> <f> + <fio/2, and for T = can be written 
as 

. , ,. , ... , 2 AhD^ . i n 

Al = I — ] — ; — > pg p sin 2np- 



P =i 



0o/2 



(44) 

In the regime of experimental interest, L 3> R 3> Ib, the 
coefficients g p , 



9p 



1 

2p 3 7T 2 



oc 
m=0 



-phj Ib yAo(«m) 



1 



IB 



(45) 

can be estimated by evaluating the sum in the saddle- 
point approximation, 



g p w 0.13p 



-3.5 



R 



^fi~B~L 



l+p 



L 

In 



All harmonics of the average persistent current are ex- 
ponentially suppressed; the higher the harmonic p, the 
stronger is the suppression. Note that this implies that 
for sufficiently strong magnetic field, measurements of the 
average current, e.g., by employing large arrays of rings, 
should be dominated by the canonical-ensemble contri- 
bution of the free-electron model 



COMPARISON TO EXPERIMENT 



The recent development of cantilever-based torsional 
magnetometers with integrated mesoscopic rings^ re- 



< o.i 



o.oi 



o.ooi 




FIG. 7: (Color online) Temperature dependence of the typ- 
ical current contribution from the p th harmonic -J (Ip(T)}. 
The markers are the data first presented in Ref. |4] The s olid 
curves represent new fits to the data using Eqs. ( 14 1 and | |40[ ) 
while the dashed curves show the fits from Ref. |4| The sam- 
ple parameters and best-fit parameters are given in Table [I] 
Closed and open markers denote measurements taken during 
different cooldowns, over different field ranges, and at differ- 
ent magnetic field orientations. In the case of the p — 1 data 
from Sample #1, the two different field ranges over which the 
closed and open markers were taken lead to slightly differ- 
ent values of the fitting function at high temperature, with 
the lower curve corresponding to the closed markers and the 
upper curve to the open markers. The new fit curves for Sam- 
ples #2 & #3 are indistinguishable, as are the old and new 
fit curves for Sample #4. 



suited in measurements of the rings' persistent current 
in the presence of large magnetic fields P Here we briefly 
review these measurements and compare them with the 
calculations from the preceding sections. This compar- 
ison is most readily performed by fitting the measured 
temperature dependence of the current to the form pre- 



dicted in Eqs. (14 1 and (40) 



The parameters characterizing each sample are col- 
lected in Table [j] T he temperature dependence of the 

p th harmonic J {I 2 ) of each sample's typical current was 

determined as follows. At a single temperature T , the 
mean square amplitude of the p th harmonic of the current 
was extracted from a measurement of 1(B) taken over a 
range of B spanning many B c . This large span ensured 
that the mean was determined from a large number of 
independent measurements, as discussed at the end of 



Sec. IIIB For each sample, the form of 1(B) was found 



to be independent of temperature except for an overall 
scaling. This scaling was determined by measuring 1(B) 
over a smaller field range (with bounds denoted by -B m i n 
and -B max ) at each subsequent temperature and compar- 
ing the magnitude of each harmonic with the value mea- 
sured over the same field range at T . This procedure, as 
well as other details of the measurements, are described 

in detail in Ref. 01 The resulting values of J (I 2 ) are 
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TABLE I: Sample parameters. "Marker" refers to the markers used in Fig. [7] with closed and open markers representing two 
different cooldowns of the same sample. For the closed markers the angle between the magnetic field and the plane of the rings 
was 6° and To = 323 mK, while for the open markers the angle was 45° and To = 365 mK. N denotes the number of rings in 
the sample. The ring circumference and linewidth are given by L and w. The thickness of each sample was 90 nm. The spin 
orbit scattering length L so = 1.1 ± 0.25 /im. B m in and -B max give the bounds for measurements of 1(B) taken over smaller field 
ranges. Dl and Dzso are extracted from fitting the persistent current data. Dl is the best-fit value of the diffusion constant 
found in Ref. 0] which assumed the limit of strong spin-orbit scattering and large Zeeman splitting. Dzso is the best-fit value 
of the diffusion constant found by taking into account the finite spin-orbit scattering rate and Zeeman splitting as described in 
Section [V] The estimated uncertainty in all fit coefficients is 6%. 



Sample 


Marker 


p 


N 


L (/»m) 


w (nm) 


-Bmin (T) 


-Bmax (T) 


D L (cm 2 /s) 


Dzso (cm 2 /s) 




• 


1 








6.2 


6.8 






#1 


O 

O 


1 

2 


1680 


1.9 


115 


5.0 
5.0 


5.2 
5.2 


271 


234 


#2 


▼ 


1 


990 


2.6 


85 


7.15 
5.39 


7.60 
5.48 


214 


195 


#3 




1 


1 


2.6 


85 


8.32 


8.40 


215 


195 


#4 


■ 


1 


242 


5.0 


85 


7.1 


7.3 


205 


196 



shown in Fig. [7j 

In Rcf. 4 , this data was analyzed by assuming the limit 
of strong spin orbit scattering: l/r so 3> {E C ,T} 7 and 
large Zeeman splitting, E z ~> {E C ,T}. As can be seen 
from the sample parameters listed in Table |TJ this as- 
sumption is fairly accurate though not exact. For these 
samples 0.075 < 1/E c t so < 0.47 while 0.15 < T/E c < 
1.7. From Fig. |6]it is clear that these parameters are not 
fully within the strong spin orbit scattering limit. Addi- 
tionally, for the smallest rings (Sample #1) the limit of 
large Zeeman splitting does not hold at the highest tem- 
peratures in Fig. [7] where T/E z ~ 0.36 and deviations 

from the large Zeeman splitting limit change by as 

much as 5%. For Samples #2, #3, and #4 T/E z < 0.17 

resulting in < 1% deviations of J (Ip) from the large Zee- 



man splitting limit. As a result the measurements of Ref. 
[|] were not fully in the strong spin orbit scattering, large 
Zeeman splitting limit, so we reanalyze the data here, 

taking into account the full dependence of y (Ip) on spin 

orbit scattering and Zeeman splitting. 

We fit the data from Ref. 0] (Fig. [7]) using the expres- 
sion for J (I 2 (T, D, L so , E z )) derived from Eqs. ( 14) and 



(40). The only fitting parameter is the electron diffu- 
sion constant D. The spin orbit length L so = yJDr so = 
l.l±0.25/im was determined independently from magne- 
totransport measurements of a wire codeposited with the 
ringsP Since each data point in Fig. [7] is extracted from 
measurements of 1(B) made over a range of B, we can- 
not use a single value of the Zeeman splitting; instead, 
we average over the magnetic field range to obtain the 
fitting function 



(Jp (-^""j -D, E SO: B m [ n: D max )) 
/ dB(P p (T, D, L SOl E Z (B))) 

V ^max — -^min 



(47) 



The best-fit values of the diffusion constant -Dzso are 
given in Table [I] The corresponding fits are shown in 
Fig. [7] as solid lines. For comparison the values of the 
diffusion constant found in Ref . [JJ D^, are also given in 
Table |T] and the corresponding fits are shown as dashed 
lines in Fig. [7J 

Figure [7] and Table [I] show that the finiteness of the 
spin-orbit scattering rate and the Zeeman energy result in 
small but noticeable changes to the fitted curves and the 
extracted values of D. We find that most of the difference 
is due to the finite spin orbit scattering rate, which leads 
to a non-negligible contribution to the current from the 
second F p term in Eq. ( 40 1 . 

The finite Zeeman energy modifies the current via the 
last two F p terms in Eq. (40), leading to a correction 
which becomes appreciable (> 1%) only for the higher 
temperature measurements of Sample #1. The resulting 
correction oscillates as a function of temperature, result- 
ing in a best-fit value of D indistinguishable from the 
case of large Zeeman splitting. 

The values of -Dzso for Samples #2, #3, and #4 
agree with each other to within the experimental uncer- 
tainty (which is estimated to be 6% in the Supplemental 
Online Material of Ref. .4). This agreement is consis- 
tent with the fact that the rings in these three samples 
have the same cross-sectional dimensions. The value of 
-Dzso = 234 cm 2 /s measured for Sample #1 is some- 
what larger, which may reflect these rings' larger cross 
section. Resistivity measurements of the codeposited 
wire having the same cross section as Sample #1 give 
D = 260 ± 12 cm 2 /s, consistent with the value measured 
for Sample #lP 



VI. CONCLUSIONS 

Motivated by a new and highly sensitive experimental 
technique^ for measuring mesoscopic persistent currents, 
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we presented a theory of persistent currents in large, but 
non-quantizing, magnetic fields. The theoretical results 
of this paper formed the basis for establishing the remark- 
able quantitative agreement between experiment and the- 
ory found in Ref. and further refined in Sec|V] To reach 
this agreement, we not only needed to take into account 
the large magnetic field, both for the single-particle and 
the interaction contributions to the persistent current, 
but also spin effects. 

In addition to forming the basis for a quantitative com- 
parison with experiment, it is also worth emphasizing 
several theoretical conclusions from our results. 

(i) The magnetic field penetrating the ring leads to 
qualitative changes in the dependence of the persistent 
current on the Aharonov-Bohm flux. At zero magnetic 
field, the persistent current is a periodic function of flux. 
Zero flux as well as integer and half-integer multiples of 
the flux quantum are special points where the persistent 
current vanishes. At large magnetic fields, the persistent 
current I(4>) is still a periodic function of flux, but the 
typical magnitude (7 2 (0)) is no longer dependent on flux. 

(ii) Previous theoretical works have shown that there 
are two principal contributions to mesoscopic persistent 
currents: a free-electron contribution and an interaction 
contribution. In experiments, it is not always easy to 
disentangle these two contributions (especially for the 
even harmonics of the persistent current). In fact, while 
the interaction contribution is expected to dominate the 
ensemble- averaged persistent current, both of them con- 
tribute significantly in single- or few-ring experiments. 
We conclude from our results that the application of a 
large magnetic field penetrating the ring strongly sup- 
presses the interaction contribution to the persistent cur- 
rent so that the technique of Ref. H] provides direct access 
to the free-electron contribution. 

(iii) One of the principal advantages of the experimen- 
tal technique of Ref. 0] is that unlike SQUID-based ap- 
proaches, it allows for measurements over a wide range 
of magnetic fields and thus of many oscillations of the 
persistent current with flux. Our results for the auto- 
correlation function of the persistent current at different 
magnetic fields imply that averaging over magnetic field 
is equivalent to an ensemble average (ergodic hypothe- 
sis). One of the possibilities raised by this result is a 
direct measurement of the entire distribution function of 
the persistent current. 

The experimental technique of Ref. |31 has brought 
many additional experiments on persistent currents and 
related phenomena within experimental reach. Our ap- 
proach should be a valuable starting point for analyzing 
such future experiments. 
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Appendix A: Arbitrary magnetic-field configurations 

Within the toroidal-field model discussed and em- 
ployed in Sees. |IIIB| and [iVj we could account for the 
magnetic field penetrating the ring by perturbation the- 
ory. This perturbative calculation was valid as long as 
R <C Ib- The toroidal-field model was special in that at 
the surface of the ring, the vector potential Ay (r) associ- 
ated with the magnetic field penetrating the ring points 
parallel to the surface. As a result, the vector potential 



does not enter into the boundary condition Eq. ( 18 ) for 
the equation of the cooperon or the diffuson. Then, com- 
puting the perturbative shift of the eigenvectors by the 
in-plane magnetic field amounts to conventional pertur- 
bation theory as familiar from quantum mechanics. 

This is no longer the case for more general (and more 
realistic) models of the in-plane field. Instead, the mag- 
netic field enters not only the diffuson or cooperon equa- 
tion, but also the boundary condition. In this appendix, 
we show how one can in principle reduce the resulting 
generalized problem of perturbation theory to the con- 
ventional case of quantum-mechanical perturbation the- 
ory. 

The equation for the cooperon or the diffuson is given 

by 



D [V - ieA ± - ieAy] 2 ip = Etp, 



(Al) 



with the appropriate choice of magnetic field. This equa- 
tion needs to be solved in conjunction with the boundary 
condition 



n ■ [V — ieAj_ — ieAn]^! 







(A2) 



valid at the surface £ of the ring. We make the gauge 
choice V • A|| =0. The basic observation is that we 
can eliminate the vector potential A|| from the boundary 
condition by the gauge transformation 



^(r) = e l/(r) ?Ai(r). 



(A3) 



The new function ifji(r) satisfies the modified diffusion 
equation 



iVf] ih=Efh. 
iV/]^| s = 0. 



- D [V — ieAj_ — zeAy 
with boundary condition 

n • [V — ieAj_ — ieAii 
If we choose the gauge transformation such that 

en - A||| s = ft- V/|jj 
combined with the gauge choice 

v 2 / = o, 



(A4) 
(A5) 
(A6) 
(A7) 
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we reduce the problem to a form which is amenable to 
standard techniques of perturbation theory, namely Eq. 



( A4 ) combined with the boundary condition 
n- [V-ieAjVils =0. 



(A8) 



The principal technical difficulty consists in solving the 
"electrostatics" problem defined by Eqs. (A6| and (A7) 
to find the function /(r). 
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